The Quark-Mass Dependence of Tc in QCD: Working up from m = or down from 

m = oo ? 



Adrian Dumitru0 Dirk R6der0 and Jorg Rupper10 

Institut fiir Theoretische Physik, J. W. Goethe Universitdt 
Postfach 11 19 32, D-60054 Frankfurt am Main, Germany 



We analyze the dependence of the QCD transition temperature on the quarlc (or pion) mass. We 
find that a linear sigma model, which links the transition to chiral symmetry restoration, predicts a 
much stronger dependence of Tc on m,r than seen in present lattice data for nin ^ 0.4 GeV. On the 
other hand, working down from m.^ = oo, an effective Lagrangian for the Polyakov loop requires 
only small explicit symmetry breaking, 61 ~ exp(— m,r), to describe rc(m^) in the above mass range. 
. . . Physically, this is a consequence of the flat potential (large correlation length) for the Polyakov loop 

I in the three-color pure gauge theory at Tc- We quantitatively estimate the end point of the line 

■ of first order deconfining phase transitions: ~ 4.2 ^/a ~ 1.8 GeV and Tc — 240 MeV for three 

' flavors and three colors. 
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I. INTRODUCTION 
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' Lattice QCD calculations at finite temperature and with dynamical fermions are presently performed for quark 
masses exceeding their physical values; for a recent review see 0. To date, pion masses as low as « 400 MeV are 
feasible 2], about three times the physical pion mass. When comparing effective theories to first- principles numerical 
data obtained on the lattice it is therefore important to fix the parameters (coupling constants, vacuum expectation 
^""^ ' values and so on) such as to match the values of physical observables, e.g. of m^, to those of the lattice calculations. 
\ For example, the QCD equation of state in the confined phase appears to be described reasonably well by that of a 
hadron resonance gas model, after extrapolating the physical spectrum of hadrons and resonances to that from the 
' lattice |3- Thus, lattice data on the dependence of various observables on the quark (or pion) mass constrain effective 
theories for the QCD phase transition at finite temperature and could provide relevant information on the driving 
O ^' degrees of freedom. 

D ', In this paper, we analyze the dependence of the chiral symmetry restoration temperature on the vacuum mass of 
• the pion using a linear sigma model in sectio njiTll The linear sigma model provides an effective Lagrangian approach 
^ I to low-energy QCD near the chiral limit 0, |Q. It incorporates the global flavor symmetry, assuming that "color" 
can be integrated out. For exaniple, it allows one to discuss the order of the Nf = 2 + 1 chiral phase transition as a 
function of the quark masses 0, la, IE • 

Instead of working up from zero quark mass, one could start with the quark masses taken to infinity, that is, with 
a pure gauge theory. Then, one can discuss the deconfinement transition at finite temperature within an effective 
Lagrangian for the Polyakov loop with global Z{Nc) symmetry H H E IH E El El IH (A^c is the number of 
colors). For finite pion mass, the symmetry is broken explicitly, and the phase transition (or cross over) temperature 
is shifted, relative to the pure gauge theory where pions are infinitely heavy. In section llVl we determine the cndpoint 
of the line of first-order transitions for three colors, and extract the magnitude of the explicit Z{3) breaking from 
lattice data on ATr. 
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II. PION MASS AND DECAY CONSTANT IN VACUUM 



The Lagrangian of QCD with the quark mass matrices set to zero is invariant under independent rotations of 
the Nf right-handed and left-handed quark fields. It exhibits a global SU{Nf)ji x SU{Nf)L symmetry, leading to 
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2{Nj — 1) conserved currents. 



Those are iV| 



1 vector currents, V^f" = qj^Xiq /2, and 



= Ql^lbKq /2, with Ai the generators of SU{Nf), normahzed according to tr A^Aj 



25^ 



1 axial currents, 
. The SUiNf)v 



subgroup of vector transformations is preserved in the vacuum [1^, while the SU{Nf)A is broken spontaneously by 
a non-vanishing chiral condensate {quqi) 7^ 0, leading to non-conservation of the axial currents. 

In reality, of course, even SU{Nf)v is broken explicitly by the non-vanishing quark mass matrix. Nevertheless, 
since at least and are very small in the physical limit, the SU{2)y symmetry is almost exact in QCD. The 
small explicit breaking of SU{2)v is responsible for the non-vanishing pion mass, as given by the Gell-Mann, Oakes, 
Renner relation 
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We neglect isospin breaking effects here, and so assume that r7i„ — rud = rriq. (qq) denotes the sum of the vacuum 
expectation values of the operators ujiUl and dadL, and their complex conjugates. The proportionality constant 
is the pion decay constant. It should be noted that is only valid at tree level, and that loop effects induce an 
implicit dependence of both /^r and (qq) on m^. For small mg, this dependence can be computed in chiral perturbation 
theory For example, at next-to-leading order. 
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where M and F are the couplings of the effective theory (equivalent to {qq) and 
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, and A3 and A4 are two 

renormalization-group invariant scales. These relations link the behavior of f-^ to that of m^, the mass of a physical 
state. (In what follows, we use TOtt to vary the strength of explicit symmetry breaking rather than using directly the 
scale dependent quark masses). 

More accurate results than eqs. H2I3|I can perhaps be obtained by computing quark propagators for various quark 
masses on the lattice. Ref. jig analyzed the propagators for gauge field configurations generated with the stan- 
dard Wilson gauge action ("quenched QCD"), using overlap fermions with exact chiral symmetry. They obtained 
a parametrisation of both tjIt^ and f-^ in terms of the mass ruq of u and d quarks (see section 2 in 18]) which 
allows us to express /^r as a function of m^. Their data covers an interval of 0.4 GeV < < 1 GcV, and 
0.15 GeV < \/2/^ < 0.22 GeV. 



III. LINEAR SIGMA MODEL AT FINITE TEMPERATURE 



In this section, we discuss chiral symmetry restoration at finite temperature, and in particular the dependence of 
the symmetry restoration temperature on the pion mass. For simplicity, we restrict ourselves here to the two-flavor 
case. Our emphasis is not on the order of the transition as the strange quark mass is varied but rather on how the 
temperature at which the transition occurs (be it either a true phase transition or just a cross over) depends on the 
pion mass. Such dependence arises from two effects. First, of course, due to explicit symmetry breaking occuring 
when TOtt > 0. Second, due to the "indirect" dependence of spontaneous symmetry breaking, i.e. of the condensate 
{qq) resp. /^, on the pion mass (through pion loops, see previous section). The tree- level potential of the linear sigma 
model with SU{2)v x SU{2)a = 0(4) symmetry is given by 

V{a,TT) = ^m^cf>' + ^^*-Ha , (4) 

with (f>a — {(J, 7?). For < the 0(4) symmetry of the vacuum state is broken spontaneously to 0(3), leading to a 
non- vanishing scalar condensate (cr) = The explicit symmetry breaking term ~ H provides a mass to the pions. 
At tree level, the masses are given by 

= rn^ + 3Acr^ , = rri^ + . (5) 

Below, we employ the Hartree approximation to investigate the dependence of the transition temperature on the 
pion mass. This approximation scheme is known to exhibit problems in the chiral limit in that the Goldstone theorem 
is violated and the phase transition is incorrectly predicted to be of first order (for Nf = 2). However, here we are 
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interested only in the model with explicit symmetry breaking, where the theory exhibits a cross over. Specifically, we 
consider the region of pion masses covered by the lattice data in 0, <; 0.4 GeV. 

In |19| | it was shown that such truncated non-perturbative resummation schemes can be renormalized with local 
counter terms obtained in the vacuum (see also "2^ for \<j)'^ theory). These ideas were applied in "21] to theories 
with global symmetries, and a BPHZ-likc renormalization scheme was introduced for the 0(4) linear sigma model in 
Hartree approximation without explicit symmetry breaking. The scheme can be straightforwardly extended to the 
case iJ > 0, see eqs. (I9lllf) below. Those renormalized gap equations coincide with those introduced first by Lenaghan 
and Rischke in ref. 22]. 

In this renormalization scheme a mass renormalization scale /i is introduced and the couplings then depend on that 
scale (cf. e.g. |23|). However, choosing 
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(6) 



the four-point coupling X{fi) = Atrcc retains its tree-level (classical) value [22 |. In other words, this renormalization 
prescription evolves the renormalization scale in such a way as to keep A constant. 
Explicitly, this leads to the following expressions for the couplings |2^ : 



^^lm|_m| ^ H = U{ml-2\fl) , ^ -± (m^ - Sm^) _ 6AQ^(m.) 
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where 
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These equations determine the couplings in vacuum in terms of TOtt, /tt and ma. The dependence of f.^ on TOtt is 
taken from the data of ref. (cf. their figs. 1, 2 and the corresponding fits therein), as mentioned above. Roughly, 
for : 0.4 GeV 1 GeV, /,r increases by about 50 %, leading to an increase of the explicit symmetry breaking 
term iJ by a factor of 10. We also require the dependence of m^ on m^, which we take from a recent computation 
with standard Wilson fermions |23. Those authors find that ma is essentially a linear function of m^. We checked 
how our results in Fig. ^ depend on this assumption by using, alternatively, a linear dependence ma — m^^ + const., 
with ma = 0.6 GeV for = 0.14 GeV. We found essentially the same dependence of Tc on m^r. 

At finite temperature, we use the effective potential for composite operators 24] to determine the masses and the 
scalar condensate in the Hartree approximation. Wc follow the derivation outlined in j2j |2^ . The resulting gap 
equations are 



H = {a)[ml-2\{a)^] 



ml = m^ + 2,\{{a)^ + [QT{ma) + Q^,{ma)] + [QT{m^) + Q^.{m^)]] 



ml = m^ + \{{af + [QT{ma) + Q^{ma)] + h[QT{m^) + Qtj.[m^)]] 
where the finite temperature contribution of the tadpole-diagram is given by 

(Pk 1 1 
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Here, ma, m^, and (tr) denote the effective masses and the scalar condensate at finite temperature, respectively. 
The self-consistent solution of the above gap equations for a given vacuum pion mass determines the temperature 
dependence of the scalar condensate as the order parameter of chiral symmetry restoration. For explicitly broken 
chiral symmetry, H > 0, the transition in this approach is a cross over. We define the cross over temperature Tc by 
the peak of d{a) /dT. 

The dependence of Tc on is depicted in Fig. ^ (left), where we have also shown lattice results obtained with 
two and three degenerate quark fiavors, respectively 2] (the Nf = 2 data with standard action, the Nf = 3 data with 
improved p4- action). Driven by the increase of both /^r and H with m^, the linear sigma model predicts a rather rapid 
rise of Tc with the pion mass, as compared to the data which is nearly flat on the scale of the figure. While lattice data 
indicate a rather weak dependence of Tc on the quark mass (see also ref. [2^1. models with spontaneous symmetry 
breaking in the vacuum naturally predict a rather steep rise of Tc with the VEV (cr) vac = fir, which itself increases 
with the quark (or pion) mass. Our findings here are in qualitative agreement with those from ref. j27| who employed 
nonperturbative flow equations to compute the effective potential for two flavor QCD within the linear sigma model. 




FIG. 1: Left: The cross over temperature Tc as a function of the (vacuum) pion mass as obtained from the hnear sigma model 
with 0(4) symmetry in comparison to lattice data for two and three flavors. The scale for both Tc and is set by the 
zero-temperature string tension in the pure gauge theory, ^Ja ~ 0.425 GeV. 
Right: The scalar condensate, (cr), as a function of temperature for various pion masses. 



They also find a steeper slope of Tc(m7r) than indicated by the lattice, even though their analysis appears to predict 
a somewhat weaker increase of /tt with ttItt than the data of [l^ . which we employ here. 

Fig. n also shows the temperature dependence of the tr condensate (right). With mg. a linear function of rn^ [23j . 
the width of the cross over is approximately independent of the pion mass for 0.4 GeV < ttItt ^ 1 GeV, while we found 
considerable broadening when mo- is linear in ttItt (not shown). The chiral susceptibility d{a)/dT at its maximum 
is « 0.25, i.e. the cross over is in fact quite broad for the range of m.^ considered. Since this is at variance with 
lattice data on QCD thermodynamics (pressure and energy density as functions of temperature, see e.g. the review 
in 0), one might argue that the cross over is in fact not driven by the order parameter field but by heavier degrees of 
freedom |17*|. Such degrees of freedom could reduce the pion-mass dependence of the transition substantially: using 
three-loop chiral perturbation theory (i.e. the non- linear model), Gerber and Leutwyler find ,17j that Tc increases 
rapidly from w 190 MeV in the chiral limit (using their set of couplings) to « 240 MeV for physical pion mass. 
However, when heavy states are included (in the dilute gas approximation), then Tc increases less rapidly, from 
K, 170 MeV in the chiral limit to w 190 MeV for physical pion mass. 

Hence, perhaps the transition is not primarily driven by an order parameter field and infrared dynamics. Another 
possible approach is discussed in the next section. 



IV. EFFECTIVE LAGRANGIAN FOR THE POLYAKOV LOOP AT FINITE TEMPERATURE 

In the limit rrig 00 the quarks decouple and drop out of the theory. The (gauge-invariant) order parameter for 
the deconfining phase transition in such a pure gauge theory with Nc colors is the Polyakov loop 

£^ ^UVe^p(^tg J^^ Ao{x,T)dT^ . (13) 

Aq denotes the temporal component of the gauge field in the fundamental representation, g is the gauge coupling, 
and path ordering is with respect to imaginary time r. Its expectation value, io{T), vanishes when T < Tc, and is 
nonzero above T^. Indeed, by asymptotic freedom, £q 1 as T ^ 00. The simplest guess for a potential for the 
Polyakov loop is: 



V{i)^-^-§\i\' + l{m' {Nc^2). (14) 
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The Polyakov loop model |i[Ti[Tl[Tl is a mean field theory for £. In a mean field analysis all coupling constants 
are taken as constant with temperature, except for the mass term, ~ About the transition, condensation of i 

is driven by changing the sign of the two-point coupling: 62 > above Tc (62(T) 1 for T ^ 00), and < below Tc- 
For two colors, H14|) is a mean field theory for a second order deconfining transition "2^. The £ field is real, and 
so the potential defines a mass: {mg/TY — [1/ Zs)d'^V / , where Zg is the wave function normalization constant for 
£ |29| |. The mass is measured from the two point function of Polyakov loops in coordinate space, oc {1/r) exp(— TO£ r) 
as r ^ 00. 

For three colors, is a complex valued field, and a term cubic in £ appears in V{t}, 

V{£) = -||£P - l^-t^ + \ (Nc = 3) . (15) 

At very high temperature, the favored vacuum is perturbative, with ~ 1, times Z(3) rotations. We then choose 
63 > so that in the Z{3) model, there is always one vacuum with a real, positive expectation value for Iq. This 
produces a first order deconfining transition, where io jumps from at T~ to ic = 2&3/3 at T+ 0,^; Tc is given by 
i>2{Tc) = ~26|/9. The i field has two masses, from its real (m^) and imaginary {fhi) parts. At T+, ^fZlmulT = l^- 
The mass for the imaginary part of i is \fZlmg{T)lT oc \/b^£; at T+, fhi/me, = 3, twice the value expected from a 
perturbative analysis of the loop-loop correlation function, obtained by expanding £ from eq. H13|l to order Aq |3(M 
This mass ratio receives corrections if five-point and six-point couplings are included in the effective Lagrangian [33 
but those are not crucial for the present discussion. We note that, in principle, all of the above coupling constants could 
be determined on the lattice. The lattice regularization requires non-perturbative renormalization of the Polyakov 
loop in order to define the proper continuum limit of £ [3ll l3^ . 

Within the above mean-field theory, dynamical quarks act like a "background magnetic field" which breaks the 
Z(3) symmetry explicitly, and a term linear in £ also appears in V{£) {33. .34 . ■35-.36j : 

V{£) = -h,^-±f- - y " y = 3, m. < 00) . (16) 

Hence, as niTr decreases from infinity, bi{m-^) turns on. The normalization of 62(2^) for T ^ cx) is such that £0 
i.e. b2{T = 00) = 1 — 61 — 63. 

We first consider the case where 61 is very small, and take the term linear in as a perturbation; then the weakly 
first-order phase transition of the pure gauge theory persists (in what follows, the critical temperature in the pure 
gauge theory with 61 = will be denoted T*). The critical temperature is determined from 



b2{Tc) = --bi^l + -^^]+0{bi) . (17) 
The order parameter jumps at Tc, from 



£oiT-) = ~ + 0{bi) , (18) 



to 



4(r/) = ^63 -^1 + 0(6?). (19) 



Note that numerically £q{T~) could be much larger than 61 if the phase transition in the pure gauge theory is weak 
and so the correlation length f — l/rrii near Tc is large (i.e. if 63 is small), as indeed appears to be the case for iVc = 3 
colors [33. In other words, it could be that on the lattice £ quickly develops a non- vanishing expectation value at 
T~ already for rather large quark (or pion) masses, but this does not automatically imply a large explicit symmetry 
breaking (see also Fig. |21 below) . 

From eq. H17|l we can estimate the shift of Tc induced by letting m-^ < 00. Writing the argument of 62 in that 
equation as T* -I- ATc and expanding to first order in ATc we obtain 

ATc 61 / 962 ,^,^.2^ 2 _ / 962 



The shift in Tc with decreasing pion mass is proportional to the expectation value of the Polyakov loop just below 
Tc] all other factors on the right-hand side of eq. H2U|) do not depend on 61 or TOtt- Numerical values for 63 and for 
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b2{T) were obtained in [Tol.[T^ l3^ by fitting the effective potential H15|) to the pressure and energy density of the pure 
gauge theory with three colors; those are 63 « 0.9 and b\/ {T*db2{T*) / dT) « 1, to within 10%. We therefore expect 
that numerically ATc/T* is roughly equal to £o{T^). 

Eqs. 1)181191) seem to indicate that the discontinuity of at Tc vanishes, i.e. that the phase transition turns into 
a cross over, at a pion mass such that 61 (m^) — 2b\/27. However, we can not really extend our 0{bi) estimates to 
the end point of the line of first-order transitions because it applies, near Tc, only if ~Ab2{Tc) ^ which translates 
into bi <Si 63/IO8, see eq. (|17ll . To find the endpoint of the line of first-order transitions we solve numerically for the 
global minimum of H16|l as a function of 62, for given &i, see Fig. |21 (left). The numerical solution is "exact" and does 
not assume that 6 1 is small. We employ &3 = 0.9 to properly account for the small latent heat of the pure gauge 
theory [lO. IIB 112) Hfl IHl • Also, for hi = Q, this 63 corresponds to £c = 0^ which is close to the expectation value of 
the renormalized (fundamental) loop for the Nc = 3 pure gauge theory js^, Ull • 

Clearly, for very small bi the order parameter £q jumps at some &2 = ^2(7c)j i-e. the first-order phase transition 
persists. (The abscissa is normalized by |62(r*)| = 2&§/9.) We find that the discontinuity vanishes at bi = 0.026(1), 
so there is no true phase transition for 61 > b^. Nevertheless, we define 62 even in the cross over regime via the peak 
of dio{b2)/db2- The shift of 62 with increasing bi can now be converted into the shift of Tc itself by expanding about 



ATc 



.552 
dT 



(21) 



T=T* 



as already discussed above. We also note that from Fig. |21 (left), the susceptibility for the Polyakov loop at its 
maximum is d£o/db2 — 3.5, 2, 1.5 for bi = 0.06, 0.1, and 0.126, respectively. That is, the cross over is rather sharp 
for the values of bi shown in the figure. 

Explicit breaking of the Z(3) symmetry of the gauge theory has previously been studied in j^H^H^, and has been 
identified as the essential factor in determining the endpoint of deconfining phase transitions. Moreover, while the 
term ~ bi quickly washes out the transition, those studies showed that along the line of first order transitions the shift 
of Tc (or, alternatively, of the critical coupling (3c) is moderate, whic h ag rees with our findings. However, the numerical 
values for the critical "external field" at the endpoint obtained in |35l Is^ from actual Monte-Carlo simulations can 
not be compared directly to our estimate for b^ because we work here with the renormalized (continuum-limit) loop, 
not the bare loop. 
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FIG. 2: Left: The expectation value for the Polyakov loop, lo{h2{T)), for various values of the explicit symmetry breaking 
coupling, b\. All curves terminate at = 1 ^ T = 00. 

Right: bi as a function of m-n, obtained by matching to three flavor lattice data for Tc(m,r). The solid line corresponds to 
an exponential increase of 61 with decreasing m-n, see text. The broken horizontal line displays the endpoint of the line of 
first-order phase transitions in terms of fei; the intersection with the fei(m-n-) curve then gives the corresponding pion mass. 

Ref. studied finite-temperature QCD with N f — 3 flavors and various quark masses on the lattice (with improved 
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p4-action), and determined the critical (or cross over) temperature as a function of the pion mass. Using eq. H21() we 
can match our ^T^jT* to the data from to determine 6i(m^). In other words, we extract the function h\(m.T^^ 
required to match the effective Lagrangian (|16|l to Tcim-^) found on the lattice. The result is shown in Fig. |21on 
the right. (Again, the pion mass is normalized to the zero-temperature string tension in the pure gauge theory, 
Vcr ~ 0.425 GcV.) 

Evidently, the sa 33% reduction of Tc from TOtt = oo (pure gauge theory) to m^r/ ^cr w 1 requires only small explicit 
breaking of the Z{3) symmetry for the Polyakov loop £: we find that bi < 0.15 even for the smallest pion masses 
available on the lattice. This is due to the rather weak first-order phase transition of the pure gauge theory with 
Nc = 3 colors, reflected by the strong dip of the string tension in the confined phase near T~ a nd of the Polyakov 
loop screening mass in the deconfined phase near T+ [sj; cf. also the discussion in [9l ll0llTTll30| . 

Moreover, 6i(m^) appears to follow the expected behavior ~ exp(— m^). The exponential fit shown by the solid 
line corresponds to 6i(m^) = aexp(— 6 TOtt/V'-'')' with a = 0.19 and b = 0.47 . Surprisingly, by naive extrapolation 
one obtains a pretty small explicit symmetry breaking even in the chiral limit, bi « 0.2. 

Finally, the endpoint of the line of first-order transitions at b^ — 0.026 (indicated by the dashed horizontal line) 
intersects the curve 6i(m7r) at mT^j^Ja « 4.2. For heavier pions the theory exhibits a first-order deconfining phase 
transition, which then turns into a cross over for rriTr ^ 4.2 yja « 1.8 GeV. According to our estimate, the endpoint of 
the line of first order transitions occurs at /^Tc/T* w 12%, which is slightly less than a previous (qualitative) estimate 
of 26% from ref. [13 ■ 

Of course, so far our analysis is restricted to pion masses 171^,/ > 1. On the other hand, one might cross a chiral 
critical point for some pion mass rriT^/ < 1 7]. Attempting a fit with the model H16|l beyond that point would then 
lead to deviations from bi ~ exp(— m^r). 

V. DISCUSSION 

Three-color QCD exhibits a (weakly) first-order deconfining phase transition at a temperature Tc/y/<J « 0.63 in the 
limit of infinitely heavy quarks (v^cr « 0.425 GeV denotes the string tension at T = in this theory). Near Tc, the 
screening mass for the fundamental Polyakov loop £ drops substantiall y l37[. and so one might hope to capture the 
physics of the phase transition using some effective Lagrangian for £ 11 lalll 111 111 111 111 111. 
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FIG. 3; Schematic phase diagram in the temperature vs. quark mass plane |2J. C is the chiral critical point, D the deconfining 
critical point. 

For finite quark masses, a term linear in £ appears which breaks the Z{3) center-symmetry explicitly. This reduces 
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the deconfinement temperature, with ATc/T* on the order of the expectation value of the Polyakov loop at , cf. 
eq. CT - 

At some point then, the line of first-order deconfinement phase transitions ends 0, 13 113 1 see Fig. 01 We 

have provided a quantitative estimate of this point, ~ 4.2 ^/a ~ 1.8 GeV and Tc ~ 240 MeV for Nf = 3 degenerate 
flavors, by m atching our effective Lagrangian for the Polyakov loop to lattice data on Tc{mT^) Assuming that 
6i oc iV/ [ai shifts "D" to ~ 1.4 GeV for Nf ^2 and to 0.8 GeV for Nf = 1. 

Going to even smaller quark (or pion) masses leaves a cross over between the low-temperature and high-temperature 
regimes of QCD. The dependence of the cross over temperature on the pion mass appears to be well described by a 
small explicit breaking of the Z(3) center symmetry, bi ^ exp(— m^), down to niTr/y/a ~ 1, which is the smallest pion 
mass covered by the lattice data of ref. 0]- On the other hand, a linear sigma model leads to a stronger dependence 
of Tc on m^n- than seen in the data. 

In turn, in the chiral limit, and for Nf = 3 flavors, one expects a first-order chiral phase transition "4, 7, 39]. The 
linear sigma model should then be an appropriate effective Lagrangian for low-energy QCD The first-order 

chiral phase transition ends in a critical point "C" if either the mass of the strange quark or those of all three quark 
flavors increase. Given that the explicit symmetry breaking term for the Polyakov loop remains rather small when 
extrapolated to 0, that is bi 0.2, we speculate that "C" might be rather close to the chiral limit. Indeed, 

recent lattice estimates for Nf — 3 place "C" at itLt, ~ 290 MeV 4(| for standard staggered fermion action and A^t = 4 
lattices; improved (p4) actions predict values as low as ~ 67 MeV |4l|. 

Of course, the question arises why, for pion masses down to ~ 400 MeV, the QCD cross over is described rather 
naturally by a slight "perturbation" of the — oo limit, in the form of an explicit breaking of the global Z{Nc) 
symmetry for the Polyakov loop. Physically, the reason is the flatness of the potential for £ in the pure gauge 
theory at see e.g. the flgures in 10, 38], which causes the sharp drop of the screening mass for ^ near T+ |37l |. 
This is natural, given that finite-temperature expectation values of Polyakov loops at Nc = 3 are close to those at 
Ac = oo js^, where the potential at Tc becomes entirely flat [s^, E2- Hence, a rather small "tilt" of the potential 
(due to explicit symmetry breaking) quickly washes out the deconfining phase transition of the pure gauge theory, 
and causes a significant shift ATc of the cross over temperature already for small bi. If so, then for Ac oo, at the 
Gross- Witten point, the endpoint "D" should be located at 6i =0; the discontinuity for the Polyakov loop at Tc, 
which in a mean-field model for the pure gauge theory is 1/2 at Ac = oo [12,^2 then vanishes for arbitrarily small 
explicit symmetry breaking. This has previously been noted by Green and Karsch Is^l within a mean-field model. If 
confirmed by lattice Monte-Carlo studies, we might improve our understanding of the degrees of freedom driving the 
QCD cross over for pion masses above the chiral critical point "C" . 
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